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Statistical properties of the transmittance (T) and reflectance (R) of an amplifying layer with one- 
dimensional disorder are investigated analytically. Whereas the transmittance at typical realizations 
decreases exponentially with the layer thickness L just as it does in absorbing media, the average 

\Q • (T) and (R) are shown to be infinite even for finite L due to the contribution of low-probable 

resonant realizations corresponding to the non-Gaussian tail of the distribution of InT. This tail 

ON ■ differs drastically from that in the case of absorption. The physical meaning of typical and resonant 

realizations is discussed. 
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■ Localization of waves in active disordered systems with gain is a subject of considerable interest given its implications 
\ for stimulated emission from random media [1]. Very recently a number of theoretical [2] - [7] and experimental [8] 
' - [11] works dealing with this subject have been puplished. A recent surprising theoretical conclusion is that in one- 
dimensional random systems amplification suppresses transmittance just in the same way as absorption does, so that 
, localization is enhanced due to amplification. In contrast to this finding in this Letter we show analytically that all 
statistical moments (T m ) of the transmittance T of a randomly layered slab with gain are infinite for m > 1. 

Propagation of radiation often exhibits surprising, nonintuitive features even in homogeneous amplifying systems . 
Whereas the intensity of radiation in an infinite active media increases exponentially along the distance of propagation, 
\ the intensity of a plain wave transmitted through a finite layer as well as the energy of a point source inside the layer 
• exponentially decreases (due to the interference of multiple internal reflection from boundaries) with the increasing 
\Q ' distance . Consider for examle, wave field u satisfying the one-dimensional stationary Helmholtz equation 
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d 2 ,2 , 



u(z)=0, (1) 



Neglecting nonlinear effects, the amplification of the medium is introduced by a constant homogeneous macroscopic 
parameter - the imaginary part of the dielectric constant (see discussions in Ref. [3], [4], [7]) e(z) = e — ij (7 > for 
amplifying medium) for < z < L (e(z) — outside the layer). . The transmittance T and reflectance R of a layer 
O ■ are expressed through the Green function Go (z, z') of the wave equation (1) as 



T = \2ikG (0, L)\ 2 , R = \2ikG (L, L) - 1| 2 , (2) 



where < z < z' < L, and 

(e~ iKZ + pe +iKZ ) (e l < z '~ L ) + pe- l < z '- L )^ 
2in (er lKh — p 2 e lftL ) 



G {z,z') = (3) 



with reflection coefficient p = (k — k) / (k + k) , n = fc-y/1 + e — ij f» k' — i/2l a , l a = (k \ 1 . 

For an arbitrary value of kl — k\J\ + e, the denominator in Eq.(^) is typically a finite (i.e. non-zero) for all 
L G (0, 00), and when L/l a » 1 the transmittance is exponentially small independently on the sign of 7. Accordingly, 
for both absorption (7 < 0) and amplification (7 > 0), 

T cx exp {—L/l a ) . 

Not only the transmission coefficient, but Green function itself also exhibits this surprising duality. Indeed, in the 
limit z/l a 3> 1 and (L — z') [l a 3> 1 (corresponding to a thick layer with both z and z' taken far away from the 
boundaries), only contributions from waves reflected from the boundaries (proportional to the reflection coefficient p 
terms in Eq. (|^)) are important in the Eq.(^), and Green function Go (z, z') (field of a point source) takes the form: 

-(ik'+l/2l a )\z-z'\ /Jk'\z-z'\ \* 

Go T . = — - ■ (4) 

— 2ik \ 2in* I 
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Thus from the intensity (either transmitted or inside the layer) there is no way of telling what kind of media, amplifying 
or absorbing, we are dealing with. 

It was shown recently [5] , [7] that an apparent dual symmetry between amplification and absorption takes place also 
in random media. In particular, if the real part of dielectric constant e(z) = Ree(z) in Eq. (Q) is a random function 
of coordinate z, and L is much larger then the characteristic scattering length l s , the average specific logarithm 
of transmittance (L _1 lnP) is a negative number both for absorbing and amplifying random media: (L _1 lnT) = 
-(l/l a + l/l s ). 

We now examine the question: is duality this similarity universal, or it is a property just of the considered quantity 
^Z. _1 mT) only? We show in this Letter that generally speaking the statistics of the intensity of radiation in an 
amplifying medium differs drastically from that in an absorbing one. In particular, in systems with gain all moments 
(T m ) with m > 1 are infinite even for finite L ^> l s , while in the case of absorption (T) oc exp (—L/l a ), and dies to 
zero as L — ► oo. 

To explore the statistical properties of the transmittance in more detail, we introduce the function 

VT(R) = (T(L) m S(R-R(L))). (5) 
The statistical moments are expressed through this function as 



poo poo 

(T m )= / dRV£(R), (R m )= / dRR m Vl 
Jo Jo 



(R). (6) 



To obtain a differential equation for V™ (R) we use the following so-called invariant imbedding equations (see [12] for 
a useful introduction) for the complex reflection and transmission coefficients r and t : 

dr ik , . , ,o dint ik , , , 

TL =^r + - < L) { l + r)\— = - < L)(l + r). (7) 
Calculation of (R) by using Eqs. (5) - (7) and employing the Furutsu - Novikov formula (see [13]) yields 

l s -^V^(R)=C m V^(R), (8) 



Cm = M R \M (i? ~ 1)2 ~ 2/3 ~ 2m (i? ~ x) ) + m2R + m (/3 ~ x) ' 

where l s is the scattering length in a passive (7 = 0) media and [3 = l s /l a with l a = (k-f) 1 being a characteristic 
amplification length in the regular (e (z) — 0) system. In derivation of Eq. (ph we have assume that e(z) is taken to be 
a Gaussian random process with zero average and (e(z)e(z')) — 2cr5(z — z') (so that l s = 2/ak 2 ) and that the rapidly 
oscillating phase of the reflection coefficient is distributed uniformly and independently on the reflectance R = |r| 2 
and the transmittance T = \t\ 2 (random phase approximation). The initial and boundary conditions for the Eq.(@) 
are: V™_^ (R) — > 6(R — 0), and continuity of P™-^o (-^) a ^ ^ ne P om t R = 1- The same equation can be obtained 
from Fokker-Plank equation of Ref. [ 0] by multiplying it by T m and integrating over T. 

Equation (^) with /3 replaced by —f3 (which corresponds to absorption), for which < R < 1, was studied in [13]. 
It was found there that the presence of (3 < leads to spectrum quantization of the operator C m . A similar result 
occurs in the present case (/3 > 0) as well. The solution of (||) may be written in terms of the eigenstates of C m : 



V?(R) = J2e- E Z n)L / l *P™(R), (9) 



C m + E^)P™{R)=Q. (10) 



where P™ (R) satisfy the equations 



Keeping only leading terms in the limit R — > 00 in Eq. ( JlOj ) , one obtains 

— R - m ) PP™ (P) = 0, P -> 00, 
OR J 

which gives for the asymptotic form of the P™ (P) , 
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(11) 



It may be easily seen from df), (§), and ([Li]) that all moments (T m ) and (R rn ) diverge for m > 1. This results from 
the fact that the distribution of transmittance although being almost log-normal near its maximum [5] , [9] , changes 
it drastically in the tail, and decreases like T~ 2 for large T: 

V L {T)otT- 2 , T^oo, (12) 

In a similar fashion, the distribution of R behaves as 

V L (R) oc R- 2 , P^ oo. (13) 

We note that the distribution function for the reflectance R of an infinite (L — > oo) amplifying random media has 
been obtained in [4], but here we have shown that the average transmittance (T) and reflectance (R) (as well as all 
higher moments of T and R) are infinite for an arbitrary finite length L>I s ,l a . 

The moments (T m ) with < to < 1 are finite, and can be calculated explicitly using a suitable modification of 
approach of [13]. Let us introduce instead of R in ( |To|) a new variable x such that R = R+{x) = (tanhx/2) 2 for 
< R < 1, and R = P_(x) = (tanhx/2) -2 for 1 < R < oo. Introducing new functions \P™ 

= V? (P±(x)) \dR ± (x) /dx\ exp ($£ (x) jl) , (14) 



<&± (a;) = - In (sinhx) + 2r7iln f i (coshx ± 1)^ =p /3 cosh a;. 

It may be shown that satisfies a Schrodinger equation with imaginary "time" S = L/l 

d_ 

~ OS 

where 



9 *m = ft™^™ S = L/l,, (15) 



nt = --^ + v?{x), (i6) 

with the potential 

V£ (x) = t (1 - sinh" 2 x) ± /? (1 - m) coshx + ^(3 2 cosh 2 x. (17) 

The initial and boundary conditions for ([jl]) follow from that for Eq.(||). The spectrum {E n } of the operator L m is 
expressed through the spectra (to) of Hamiltonians 7i™ (T-C^ip™ = (to) i/CO as 

{E n } = {e+}u{e-}. (18) 

Note that if E n = e~ </ {e„} the corresponding function P™ (R) (see Eq.([|)) is identically equal zero for < R < 1. 
As the random phase approximation is expected to be valid for small [3 [14], proceeding along lines of [13] we obtain 
for In (1/(3) ^> 1 the equation 



2lfe ± r 2 (i + 2»fc±) r(i/2±(i-TO)-^j 
(i - 2*jfe±) r (1/2 ± (i - m) + ^ ' 



One can see from Eq.(|l9|) that for integer to, fc+ = k~. This property is valid for arbitrary also. For to = it 
results from the representation H± = X±X± with x± = d/dx+ (1/2) (d$j_/dx) . It is easy to see that X+X+ — X-X- 
(both TiP + and Tt°_ have the same supersymmetric partner), whence it follows (see [15]) that the spectra of TL j_ are 
identical except for the ground states. For to = 1 we have Ti\ = H°_. 
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For real (corresponding to > 1/4) and 1 < n <C ln(l//3) Eq.(fil|) gives k^ w 7m/ In (1//3), being independent 
on m in the leading on l/ln(l//3) approximation. For n = m = we have €q [m = 0) = with exact ground state 
< = V^e^exp (-$°/2) corresponding to stationary distribution V L ^oo(R) = W{R - 1)~ 2 exp (-2/3/ (i? - 1)), 
found in Ref. [4]). For small, but finite m Eq.(|l9|) yields 

e ( m ) = m(l + /?) +m 2 (1 - 2(3 In 1//3) + 0(m 3 ) (20) 

The term linear in m in this expansion can be easily reproduced for arbitrary j3 by perturbation theory: €q (m) = 

(ipo \ v - - v -\ V>o) = m (1 + P)- 

li L/l s ^> 1, only contribution from the lowest level E = £q in the sum in (9) is important, and from ( p0[ ) and (j^) 
we obtain for < m < 1 : 

(T m ) oc exp (-e (m)L/l s ). (21) 
Taking into account that (In™ T) = (d n (T m ) /dm n ) m ^ , and using Eqs. @ and @ we get: 

(InT) = -(1 -P)L/l s , var(lnT) = 2 (1 - 2/31n 1//3) , (22) 

in agreement with the results of [7] (see also [16]). 

Formulas ( |2^ ) correspond to the exponential decay of the transmittance of a random amplifying layer that at 
first glance seems to contradict to the infinite value of (T) obtained above. In fact these two results are actually 
in agreement. Since (in accordance with Eq.(p2|)) fluctuation of £ = L _1 lnT decays to zero when L — > oo, £ is a 
self-averaged random quantity whose value at each random realization tends (with the probability one) to the non- 
random limit £l— oo = — (1 — P)/l s which corresponds to the maximum of the distribution function P(£). Therefore 
the exponentially decreasing function e^ lnT ^ = e~( 1 ~^ L / ls represents the behavior of the random quantity T at the 
most probable (typical) realizations. At the same time the main contribution to the average transmittance (T) comes 
from rare, low-probable (so-called representative) realizations (corresponding to the far tail of the distribution function 
P {£))■> where the transmittance takes (due to resonances) an infinite value. 

The existence of two types of realizations (typical and resonant) can be demonstrated by the simplest example 
of a random ensemble of amplifying systems: Fabry-Perot resonators (dielectric layers) with homogeneous at each 
realization, but randomly fluctuating from sample to sample width L and real part of dielectric constant s(z). As 
was shown above, the transmittance of a dielectric layer is exponentially small independent of the sign of 7, i.e. both 
for absorption and amplification. However, in the case of amplification (7 < 0) there exist a discrete set of resonant 
pairs {k' n , L n } for which the denominator of Eq. (^) equals zero, and Go (z, z') as well as the transmittance T and 
reflectance R become infinite. The occurrence of such resonant realizations may give rise to the divergences of (T) 
and (R) . Mathematically, if we first average T (and R) over the rapid phase oscillations related to the wave length 
2-7r/fc' <C l a (this is equivalent to the random-phase approximation used in the derivation of Eq.(^|)), and then average 
the resulted "smoothed " functions over random e, we obtain (due to the infinite contributions of the above-mentioned 
resonant realizations) infinite values for (T) and (R). 

To conclude, statistical moments of the transmittance (T m ) and reflectance (R m ) for a disordered amplifying layer 
of the thickness L are calculated analytically. When m < 1 moments (T m ) as well as (InT) and varlnT exhibit 
the same L-dependance as in absorbing (Im e — > — Im e) media. This behavior is inherent for typical realizations 
corresponding to the Gaussian- like part of the InT distribution function near its maximum. However, the duality 
between absorbing and amplifying media breaks down for m > 1: all moments (T m ) and (R m ) with m > 1 diverge 
(due to the contribution of low-probable resonant realizations related to the non-Gaussian tail of the distribution) for 
finite L that is larger than the characteristic scattering length l s . 

We are grateful to I. Freund for stimulating discussions and useful criticism. 
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(lnT) can be derived directly from Eq. ([?]). After a simple algebra one gets l a dlnT/dL = l B d]n \R — 1| /dL — 
(3 |1 + r\ 2 I (R — 1). Neglecting oscillating terms and averaging with stationary distribution of R, one obtains (h/L) (lnT) = 
(p (R + 1) /(«-!)) = -(l + /3). 
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